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EFFECTS OF TURBULENT TRANSFER ON THE CRITICAL BEHAVIOUR 
N. V. Antonov, A. S. Kapustin, A. V. Malyshev 

Abstract. Critical behaviour of two systems, subjected to the turbulent mixing, 
is studied by means of the field theoretic renormalization group. The first system, 
described by the equilibrium model A, corresponds to relaxational dynamics of a 
non-conserved order parameter. The second one is the strongly nonequilibrium 
reaction-diffusion system, known as Gribov process or directed percolation process. 
The turbulent mixing is modelled by the stochastic Navier-Stokes equation with 
random stirring force with the correlator cx 6{t — t')p^~'^~y , where p is the wave 
number, d is the space dimension and y the arbitrary exponent. It is shown that, 
depending on the relation between y and d, the systems exhibit various types of 
critical behaviour. In addition to known regimes (original systems without mixing and 
passively advected scalar field), existence of new strongly nonequilibrium universality 
classes is established, and the corresponding critical dimensions are calculated to the 
first order of the double expansion in y and e = A — d (one- loop approximation) . 



1. Introduction 

Various systems of very different physical nature exhibit interesting singular behaviour 
in the vicinity of their critical points. Their correlation functions acquire self-similar 
form with universal critical dimensions: they depend only on few global characteristics 
of the system (like symmetry or space dimension). Quantitative description of critical 
behaviour is provided by the field theoretic renormalization group (RG). In the RG 
approach, possible types of critical regimes (universality classes) are associated with 
infrared (IR) attractive fixed points of renormalizable field theoretic models. Most 
typical equilibrium phase transitions belong to the universality class of the 0„-symmetric 
ip'^ model of an n-component scalar order parameter. Universal characteristics of the 
critical behaviour depend only on n and the space dimension d and can be calculated 
within various systematic perturbation schemes, in particular, in the form of expansions 
in e = 4 — (i or 1/n; see the monographs [H Ej and the literature cited therein. 

Aleksandr Nikolaevich Vasiliev made valuable contribution to the development of 
field theoretic methods and their application in the theory of critical behaviour and 
theory of turbulence. His work in this field is summarized in the three monographs 
[2], [31 H]. The most remarkable specific achievements are probably the calculation 
of Fisher's exponent rj in the (9„-symmetric ip'^ model to the order l/rr" [5] and the 
third-order calculation of the anomalous exponents in Kraichnan's model of turbulent 
advection [6]. In the present paper, we apply the field theoretic RG and generalized e 
expansion to the problem of the effects of turbulent transfer on various types of critical 
behaviour. 

Over the past few decades, constant interest has been attracted by the spreading 
processes and corresponding nonequilibrium phase transitions; see e.g. the review papers 
[3, [8] and the literature cited therein. Spreading processes are encountered in physical. 
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chemical, biological and ecological systems: autocatalytic reactions, percolation in 
porous media, epidemic diseases and so on. The transitions between the fluctuating 
(active) and absorbing (inactive) phases, where all the fluctuations cease entirely, are 
especially interesting as examples of nonequilibrium critical behaviour. 

It has long been realized that the behaviour of a real critical system is extremely 
sensitive to external disturbances, gravity, impurities and turbulent mixing; see the 
monograph p] for the general discussion and references. What is more, some 
disturbances (randomly distributed impurities or turbulent mixing) can produce 
completely new types of critical behaviour with rich and rather exotic properties. 

These issues become even more important for nonequilibrium phase transitions, 
because the ideal conditions of a "pure" stationary critical state can hardly be achieved 
in real chemical or biological systems, and the effects of various disturbances can never 
be completely excluded. In particular, intrinsic turbulence effects cannot be avoided 
in chemical catalytic reactions or forest flres. One can also speculate that atmospheric 
turbulence can play important role for the spreading of an infectious disease by flying 
insects or birds. Effects of different kinds of regular and turbulent flows on the critical 
behaviour were studied in [T0]-p!7]. 

In a number of papers [ll]-[I7j, critical behaviour of various systems, subjected 
to the turbulent mixing, was studied by means of the fleld theoretic RG. As a rule, 
the turbulence was modelled by the time-decorrelated Gaussian velocity fleld with 
the velocity correlation function of the form (vv) oc 5{t — t')p~'^~^, where p is the 
wave number and < ^ < 2 is a free parameter with the real ( "Kolmogorov" ) value 
^ = 4/3. This "Kraichnan's rapid-change model" has attracted enormous attention 
recently because of the insight it offers into the origin of intermittency and anomalous 
scaling in fully developed turbulence; see the review paper |18] and references therein. 
The RG approach to that problem is reviewed in pj]. In the context of our study it is 
especially important that Kraichnan's ensemble allows one to easily model anisotropy 
of the flow [fT6] and compressibility of the fluid [17] , which appears much more difficult 
if the velocity is described by the full-scale dynamical equations. 

However, the Gaussianity and vanishing correlation time are drastic simplifications 
of the real situation, and it is desirable to investigate effects of turbulent mixing, caused 
by more realistic velocity fields. In this paper, we study effects of a strongly non- 
Gaussian velocity field with finite correlation time, governed by a stochastic dynamical 
equation. More precisely, we employ the stochastic Navier-Stokes equation for an 
incompressible velocity, with random stirring force with the correlator cx p^~'^~y^ where 
y is the arbitrary exponent with the physical ("Kolmogorov") value y = 4. The RG 
approach to this model is reviewed in [21 H] . 

Two representative cases of dynamical critical behaviour are considered: 
equilibrium critical dynamics of a non-conserved order parameter with ip'^-tjpe 
Hamiltonian, and the nonequilibrium system near its transition point between the 
absorbing and fiuctuating states. The former model corresponds to critical fiuid systems 
(binary mixtures or liquid crystals), and the latter describes the spreading processes in 
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reaction-diffusion systems, belongs to the universality class known as Gribov process or 
directed percolation process, and is equivalent (up to the Wick rotation) to the Reggeon 
field theory 

It is shown that, depending on the relation between y and d, the both systems 
exhibit various types of critical behaviour, associated with different IR attractive fixed 
points of the RG equations. In addition to known asymptotic regimes (like equilibrium 
critical dynamics without mixing or passively advected scalar without self- interaction), 
existence of new, strongly nonequilibrium, types of critical behaviour (universality 
classes) is established, and the corresponding domains of stability in the y-d plane 
and the critical dimensions are calculated to the leading order of the double expansion 
in y and e = A — d, which corresponds to the one-loop approximation of the RG. 

2. Description of the models 

In the Langevin formulation the models are defined by stochastic differential equations 
for the order parameter ip = ip{t, x): 

dt^ = Ao {(-To + d')^ - uo^ysi] + C (2.1) 

for the model A and 

dt^ = Ao {(-To + d')^P - + Cv^ (2.2) 

for the Gribov process. Here, dt = d/dt, is the Laplace operator, qq and Uq > are 
the coupling constants, Aq > is the kinematic (diffusion) coefficient and tq oc (T — Tc) 
is the deviation of the temperature (or its analog) from the critical value. The Gaussian 
random noise ( = C{t,x) with zero average is specified by the pair correlation function: 

(C(t,x)C(t',x')) =2Ao5(t-t')'5^'Hx-x') (2.3) 

for the model A and 

(C(t,x)C(t',x')) =^7oAo5(t-t')5^'Hx-x') (2.4) 

for the Gribov process; d being the dimension of the x space. The factor in the 
noise term of (12. 2p guarantees that in the absorbing state the fluctuations cease entirely. 
The expressions for the correlators differ only by normalization: the factor 2Ao in (12. 3 p 
is dictated by the fluctuation-dissipation relation and ensures the correspondence to 
the static ip'^ model, while QqXq in (12. 4p provides the simple form of the symmetry that 
exists in the field theoretic formulation of the Gribov model; see eq. (13. 5 p below. The 
subscript "0" marks the bare (unrenormalized) parameters; their renormalized analogs 
(without the subscript) will appear later. 

For incompressible fluid, the Galilean covariant coupling with the transverse (due 
to the incompressibility condition diVi = 0) velocity field v = {t'j(t,x)} is introduced by 
the substitution 



dt ^Vt = dt + Vidi 



(2.5) 
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in fl2.ip and fl2.2p . where di = d/dxi and Vj is the Lagrangian (Gahlean covariant) 
derivative. We will employ the velocity field satisfying the NS equation with a random 
stirring force 

VtVk = Jyod'^Vk - dkV + fk, (2.6) 

where is the Lagrangian derivative ( 12. 5 p , V and fk are the pressure and the transverse 
random force per unit mass. We assume for / a Gaussian distribution with zero average 
and correlation function 

{M^)fA^')) = %^ [ rfpP.,(p)^^/(p)exp{ip(x-x')}, (2.7) 

l^^J Jp>m 

where Pij{p) = Sij —PiPj/p'^ is the transverse projector and 'Df{p) is some function of 
p = IpI and model parameters. The momentum m = 1/C, the reciprocal of the integral 
turbulence scale £, provides IR regularization (its precise form is unessential; the sharp 
cutoff is the simplest choice for the practical calculations). 

The standard RG formalism is applicable to the problem (12. 6p . (12. 7p if the 
correlation function of the random force is chosen in the power form 

Vf{p) = Dop'~'~\ (2.8) 

where > is the positive amplitude factor and the exponent < y < 4 plays the 
role of the RG expansion parameter, analogous to that played hj e = 4 — d in models 
of critical behaviour. Its physical value is y = 4: with the appropriate choice of the 
amplitude, the function (12. 8 p for y — )■ 4 turns to the 6 function, T>f{p) oc 5(p), which 
corresponds to the injection of energy to the system owing to interaction with the largest 
turbulent eddies; for a more detailed discussion of this point see e.g. [2111]. 



3. Field theoretic formulation and renormalization 

The stochastic problems (I2.ip -f l274l) can be reformulated as field theoretic models of the 
doubled set of fields \E' = {ipjifj"^} with action functional 



5^(V', ^^) = {-dt + \,d^ - Aoro) ^ + Ao^^t - V (3-1) 
for the model A and 



5g(^, i^^) = i^\-dt + \,d^ - Aoro)^ + ^ {(V^t)2^ _ ^t^2| (3 2) 

for the Gribov model. Here, ip'^ = '?/'^(t, x) is the auxiliary "response field" and the 
integrations over the arguments of the fields are implied, for example 



iP^dt'4) = j dt j dxip''{t,x)dti^{t,x) 



The stochastic problem (I2.6l) - (l2.8p corresponds to the field theoretic model with 
the action 

SNsiv\v) = v'Dy/2 + v' {-Vt + uod^}v, (3.3) 
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where is the correlation function f l2.7p and all the needed integrations and 
summations over the vector indices are understood. The auxiliary vector field v' = 
x)} is also transverse, div[ = 0, which allows one to omit the pressure term in the 
action functional fl3.3p . 

The field theoretic formulation means that statistical averages of random quantities 
in the original stochastic problems can be represented as functional integrals over the 
full set of fields with the weight expi5($), and can therefore be viewed as the Green 
functions of the field theoretic models with actions f l3.ip -( l3l6l) . In particular, the linear 
response function of the stochastic problems (I2.ip - fl2.4p is given by the Green function 

G= (^^(t,x)V^(t',x')) = Jvij^ j Vi) V^t(t,x)V'(t',x') exp5(V^,V^t) (3.4) 

of the corresponding field theoretic models. 

The model (13. 2 p is symmetric with respect to the transformation 

V^(t,x) ^ V^'^(-t,-x), V^'f(t,x) V^(-t,-x), ^o^-^o- (3.5) 

Refiection of the constant is in fact unimportant because the actual expansion 
parameter in the perturbation theory of the model is uq = Qq. The model (13. ip is 
symmetric with respect to the refiection of the fields ip — > —ip, ip'^ — )■ —ip'^. These 
symmetries survive the inclusion of the velocity field. 

The full-scale models are described by the action functionals 

Slam=SNs{v\v)+SA,G{^,^\v), (3.6) 

where \& = {tp^ip'^ ,v,v'} is the full set of fields and the substitution (12. 5 p is made in 
the functionals ( 13. ip and ( 13. 2p . For these models, the full set of coupling constants 
("charges") involves the three parameters 

Mo ~ A'"^ = Do/u^ ~ eo = Aq/z^o, (3.7) 

where A is some typical UV momentum scale. The ratio Cq is not an expansion parameter 
in the perturbation theory, but it should also be treated as an additional coupling 
constant because it is dimensionless and the renormalization constants and RG functions 
depend on it. 

From the relations (13.70 it follows that the interactions ('?/'^)^'?/' and ip^ip"^ in (13.20 
and ip'^'ip^ in (13. ip become logarithmic (the corresponding coupling constant uq becomes 
dimensionless) at d = 4. Thus for the single-charge problems ( 13. ip . (13.20 . the value 
d = 4 is the upper critical dimension, and the deviation e = A — d plays the part of the 
formal expansion parameter in the RG approach: the critical exponents are nontrivial 
for e > and can be calculated as series in e. The vertex term v'{vd)v in (13.30 and the 
additional interactions tp'^{vd)ilj in the full models (13. 6p become logarithmic at y = 0. 
The parameter y is not related to d and can be varied independently. However, for 
the RG analysis of the full problems it is important that all the interactions become 
logarithmic at the same time. Otherwise, one of them would be weaker than the others 
from the RG viewpoints and it would be irrelevant in the leading-order IR behaviour. 
As a result, some of the scaling regimes of the full model would be lost. 
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Figure 1. One- loop approximation for the relevant 1-irreducible Green functions in 
the model 

In order to study all possible scaling regimes and the crossovers between them, we 
need a genuine three-charge theory, in which all the interactions are treated on equal 
footing. Thus we will treat e and y as small parameters of the same order, e oc y. Instead 
of the plain e expansion in the single-charge models, the coordinates of the fixed points, 
critical dimensions and other quantities will be calculated as double expansions in the 
plane around the origin, that is, around the point in which all the coupling constants 
in (13.71) become dimensionless. 

The analysis based on the dimensionality considerations and the symmetries of 
the full-scale models (13. 6p shows that they are multiplicatively renormalizable. The 
role played by the symmetries is very important: in particular, the Galilean invariance 
requires that the counterterms ip^dtip and il!^{vd)ijj enter the renormalized action only 
in the form of invariant combination ip^Vt'ip- It also shows that the counterterm ip^ipv'^, 
absent in the unrenormalized actions (13. 6p and allowed by the dimension, is in fact 
forbidden. Thus, all the UV divergences (having the form of singularities at e and ?/ — )■ 0) 
can be absorbed into a finite set of renormalization constants Zi. The renormalized 
action functionals have the forms: 

= S§s + { -Zidt + Z^\d^ - ZsAr} ^ + Z^XiIj^^ - ^5^^ V'(3.8) 
for the model A and 

= ^_z^dt + Z^Xd^ - Z^Xt] ^P+Z,^^ {{^^fij - tAV'} (3.9) 

for the Gribov model, where S§g is obtained from (13.31) by the substitution uq vZy 
and WqV^ — )• wfi^u^ (the nonlocal term with the random force correlator in (13.31) is 
not renormalized). Here and below r, u and so on are renormalized analogs of the 
bare parameters (with the subscripts "0") and fi is the reference mass scale (additional 
arbitrary parameter of the renormalized theory). 

The one-loop calculation of the renormalization constants is easily performed: 
in fact, in this approximation there are no new Feynman diagrams in comparison 
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to the models fl3.ip - fl3.3l) and the passive scalar case. More precisely, all the new 
diagrams appear UV finite for an incompressible fiuid and give no contribution to 
Zi. This fact is illustrated for the model A in figure [H where all 1-irreducible Green 
functions, needed for the calculation of the renormalization constants, are shown in 
the one-loop approximation. The solid lines with arrows denote the bare propagator 
the arrow pointing to the field ip'^ . The solid lines without arrows correspond 
to the propagator {ipip)o and the wavy lines denote the velocity propagator {vv)o. The 
external ends with incoming arrows correspond to the fields ip"^ , the ends without arrows 
correspond to ip. The quartic vertex with one incoming arrow corresponds to the 
interaction —ufi'^Xtp'^ilj^/3\, while the triple vertex with one wavy line corresponds to 
—il)\vd)ilj. Due to the transversality of the velocity field, the derivative at the latter 
vertex can also be moved onto the field ip"^ using integration by parts: —ip^{vd)ip = 
ip{vd)%p'^ . Thus in any diagram involving n external vertices of this type, the factor 
p" with n external momenta p will be taken outside the corresponding integrals. This 
reduces the dimension of the integrand by n units and can make it UV convergent. 
In the case at hand, this proves the UV finiteness of the last two diagrams in the 
function [ip'^ ipipip) i^ir and the only diagram in ((ip'^'ip'^)i^ir (which otherwise would be 
logarithmically divergent). Since in our models the scalar field is passive (no feedback 
on the velocity statistics), the constant Z^ in the full model (13. 6p is the same as in (13. 3p . 
but with the substitution d = 4. 

In the minimal subtraction scheme the one-loop expressions for the constants Zi 
contain only simple poles in e and y and have the forms: 

Z,=Z. = 1. Z3 = l + ^. = = l + f (3.10) 



for the model A, 



11 u u 

Zi = l + -, ^3 = 1 + -, z, = i + -, 

Ae 2e e 

Z, = l + ^- 7 , (3.11) 

8e %e(e +1) ^ ' 

for the Gribov model and 

Z^ = l-— (3.12) 
I2y ^ > 

for the both cases (in order to simplify the coefficients, the factor l/lGvr^ is absorbed 
into the constants u and w). 



4. Fixed points and scaling regimes 

The RG equations for our multiplicatively renormalized models (13. 8p . (13. 9p are derived 
in a standard fashion, similar to that for the analogous models with Kraichnan's velocity 
field [T7], and we do not present them here. 

It is well known that possible IR scaling regimes of a renormalizable field theoretic 
model are associated with IR attractive fixed points of the corresponding RG equations; 



8 



see e.g. [I1I2]- For a given point, the Green functions demonstrate self-similar (scaling) 
asymptotic behaviour in the IR range, with definite critical dimensions of all fields 
and parameters F of the model. The coordinates gi* of the fixed points are found from 
the requirement that the /3-functions, corresponding to all renormalized couplings Qi, 
vanish. The type of a fixed point is determined by the matrix VLik = df3i/dgk, where (3i 
is the full set of /3-functions and Qk the full set of couplings. For an IR attractive fixed 
point the matrix Q is positive, i.e., the real parts of all its eigenvalues are positive. 

In our case, Qi = {u,w,e}. Admissible fixed point must be IR attractive for some 
values of y and e and satisfy the conditions u^:,w^,e^: > 0, which follow from the physical 
meaning of these parameters. The functions Pi, calculated in the one-loop approximation 
from the renormalization constants fl3.10p - fl3.12p . have the forms: 

for the model A, 

= + _ + (4.2) 

for the Gribov model and 

Pn^ = w {-y + w/4} (4.3) 

for the both models, with higher-order corrections in u and w. 

The analysis of the functions fl4.ip . fl4.3p reveals four admissible fixed points of the 
model A: 

(1) The Gaussian (free) fixed point: = = 0, e^, arbitrary. This point is IR 
attractive for y < 0, e < 0. The critical dimensions are found exactly: 

A^ = rf/2-l, A^t=rf/2 + l, A^ = A, = 2. 

(2) The point = 0, to* = Ay, 2e* = — 1 + a/TS (the positive root of the equation 
e(e+l) = 3), corresponding to the passively advected scalar without self-interaction: the 
vertex ip'^ip^ in (13. II) is IR irrelevant in the sense of Wilson. This point is IR attractive 
for y > 0, y > 36/2. The critical dimensions are also known exactly: 

A^ = rf/2-l, A,^t=rf/2 + l, A^ = A, = 2-y/3. 

(3) The point = 0, = e/3, arbitrary, corresponding to the pure A model: 
the turbulent advection is IR irrelevantly This point is IR attractive ioi y < 0, e > 0. 
The critical dimensions for this regime depend only on e: 

A^ = 1 - e/2, A^t = 3 - e/2, A^ = 2, A, = 2 - e/3, 

with the higher-order corrections, known up to for A^ [20] and for the others 

I This becomes obvious if, by rescaling the fields, the coupling constant wq is placed in front of the 
interaction terms ip'{vd)'ip, which is more familiar for the field theory. We do not do it, however, in 
order to retain the natural form of the covariant derivative (|2.5p . 
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(4) The most interesting point = Ay, = e/3 — 2y /9, 2e* = — 1 + y/TS, IR 
attractive for y > 0, y < 3e/2. It corresponds to a new full-scale nonequilibrium 
universality class, in which both the self-interaction and turbulent mixing are relevant. 
Here, the critical dimensions are calculated in the form of double series in e and y. The 
one-loop expressions read: 

A, = l-l^i±i^). A„=3-i<i±i^. A. = 2-^4 A„ = 2-|.(4.4) 

The last dimension is exact, the others have higher-order corrections in e and y. 

For the Gribov case, the analysis of the functions fl4.2l) . (14. 3p reveals five admissible 
fixed points: 

(1) The Gaussian point: = = 0, arbitrary, attractive for ?/ < 0, e < 0. 
Here, the critical dimensions are: 

A^, = A^t = d/2, A^ = Ar = 2. 

(2) The point = 0, = Ay, 2e* = — 1 + y/TS, attractive for y > 0, y > 3e/2. 
It corresponds to the passively advected scalar without self-interaction: the vertices 
('?/;''') ^t/; and ip'^ip'^ in (13. 2 p are irrelevant. The critical dimensions are: 

A^ = A^t = d/2, A^ = Ar = 2- y/3. 

(3) The point = 0, = 2e/3, e* = oo, IR attractive for y < 0, e > 0. It 
corresponds to the pure Gribov process (turbulent advection is irrelevant). The critical 
dimensions depend only on e: 

A^ = A^t = 2 - 75/12, A, = 2-5/4, A^ = 2-5/12, (4.5) 

with known corrections of order [8]. 

(4) The full-scale point, corresponding to a new universality class: = Ay, 
= Ae/5 — 8y/15, 2e* = — 1 + \/l + AOy/ {Ay — e). This point is IR attractive for 

y > e/A, y < 3e/2. The dimensions are calculated as double series in e and y with the 
one-loop expressions: 

A^ = A^t =2-^ + ^, A. = 2-^, A^ = 2- I (exact). (4.6) 

(5) The point u^^ = 2e/3, = Ay, = oo, IR attractive for y > e/A, 5 > 0. 
This point requires a careful interpretation. Although the value of at this point 
is nontrivial (and the velocity field is therefore non-Gaussian), the turbulent mixing is 
nevertheless irrelevant. Indeed, straightforward analysis of the Green functions with the 
scalar fields ip, ip'^ , for example the function (13. 4p . shows that the Feynman diagrams 
involving the velocity field vanish in the limit e = A/z/ — t- oo, while the diagrams 
without the velocity are independent of e and remain finite. Thus from the physics 
viewpoints this regime is similar to (3) and corresponds to the pure Gribov process, and 
the corresponding critical dimensions indeed coincide with (14. 5p . 

In figure [2] we show the domains in the e-y plane, where the fixed points listed 
above are IR attractive. The plot corresponds to the Gribov case; for the model A, the 
domain (5) is absent, and the boundary between the domains (3) and (4) is given by 



10 

y 









y = 2£/3 


o 

z 


o 

I 


/ 4 


5 


?/ — 


m 




3 




2/ = 



e/4 



Figure 2. Domains of IR stability of tlie fixed points in the models p.6p . The 
numbers in boxes correspond to the fixed points (l)-(5) in the text. For the model yl, 
the domain (5) is absent, and the boundary between (3) and (4) is given by the ray 
?/ = 0, £ > 0. 

the ray y = 0, e > 0. In the one-loop approximation, all the boundaries of the domains 
are given by straight lines; there are neither gaps nor overlaps between the domains. 
Due to higher-order corrections to the functions f l4.ip -( H73|) . the boundaries between the 
domains (2), (4) and (5) can change and become curved. It can be argued, however, that 
no gaps nor overlaps can appear between them to all orders; cf . [171 El] ■ It is important 
here that the special cases u = or w = of the full models are "closed with respect 
to renormalization" in the sense that the functions fi^ for w = coincide with the fi 
functions of the Gribov model or model while the functions fi^^e for m = coincide 
with their counterparts in the passive scalar model to all orders of the perturbation 
theory. It is also not impossible that the absence of the regime (5) for the model A is 
an artefact of the one-loop approximation, and it will appear on the two-loop level due 
to nontrivial contributions to the renormalization constants Zi 4 in fl3.10p . 

5. Conclusion 

Effects of turbulent mixing on the critical behaviour were studied. Two representative 
models of dynamical critical behaviour were considered: the model ^4, which describes 
relaxational dynamics of a non-conserved order parameter in an equilibrium critical 
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system, and the strongly nonequilibrium Gribov model, which describes spreading 
processes in a reaction-diffusion system. The turbulent mixing was modelled by 
the stochastic Navier-Stokes equation with random stirring force with the prescribed 
correlation function oc S{t — The original stochastic problems can be 

reformulated as multiplicatively renormalizable field theoretic models, which allows one 
to apply the field theoretic RG to the analysis of their IR behaviour. We showed 
that, depending on the relation between the spatial dimension d and the exponent 
y, the models exhibit different critical regimes, associated with possible IR attractive 
fixed points of the RG equations. For the both models, the most interesting point 
corresponds to a new type of critical behaviour, in which the nonlinearity and turbulent 
mixing are both relevant, and the critical dimensions depend on the two parameters d 
and y. Practical calculations of the dimensions and the domains of IR stability for all 
the regimes were performed in the one- loop approximation of the RG, which corresponds 
to the leading order of the double expansion in y and e = 4 — d. 

From the dimensions of the coupling constants (13.71) one could expect that the full- 
scale regime (4) must take place when y and e are both positive, but the careful RG 
analysis has shown that the domains of its IR stability is in fact much narrower: for the 
Gribov model, in the one- loop level it reduces to the sector e/A < y < 2e/3, while for the 
model A one obtains < y < 2e/3. This effect leads to interesting physical prediction: 
in contrary to what could be naively anticipated, the most realistic spatial dimensions 
(i = 2 or 3 and the Kolmogorov exponent ?/ = 4 for the fully developed turbulence lie in 
the domain of IR stability of the passive-scalar regime. For the Gribov case this means 
that the spreading of the agent is completely determined by the turbulent transfer. 
For the equilibrium model A, this is reminiscent of the observation made in [TOl [TT] 
(however, for a conserved order parameter and non-random velocities) that the critical 
fluctuations are suppressed by the motion of the fluid and the behaviour of the system 
becomes close to the mean-field limit in a strong shear flow; see also discussion in [13]. 

It is interesting to compare our results with those, obtained earlier in [17], where 
the turbulence was modelled by Kraichnan's ensemble - the time-decorrelated Gaussian 
velocity field with the correlator oc 6{t — t') p~'^~^. It turns out, that the number and the 
character of the critical regimes (free theory, passive scalar, ordinary phase transition 
and the new full-scale regime) are the same for the both ensembles. (For the Gribov 
case, the single passive-scalar regime for Kraichnan's ensemble corresponds to the set 
of two regimes (3) and (5) in the Navier-Stokes model.) What is more, in the one-loop 
approximations the domains of IR stability in the y-e plane and the explicit expressions 
for the critical dimensions coincide for the two ensembles. (To compare the results for 
the two different ensembles, one has to identify y = 3^, because ^ = 4/3 for Kraichnan's 
ensemble and y = 4 for the Navier-Stokes case correspond to Kolmogorov's velocity 
spectrum.) Such agreement allows one to conclude that Kraichnan's ensemble, in spite 
of its relative simplicity, may serve as acceptable model of turbulent mixing. 
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